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ABSTRACT 


Q-mode cluster analysis of samples of the living fauna and 
of shelly remains reveal groups of samples that correspond to 
the major habitats'of Bahia la Choya. Analysis of samples ol 
the live fauna defines four groups: inner flats and marsh, outer 
and middle sand flats, rocky flats and channel, and subtidal. 
Analysis of samples of the shelly remains defines five groups: 
inner flats, marsh and channel, sandy subtidal, outer and 
middle sand flats, rocky subtidal, and rocky flats. Post-mor¬ 
tem transportation causes the association of channel samples 
with marsh and inner flat samples of shelly remains. Extensive 
time-averaging and the accumulation of diagnostic species in 
the more diverse samples ol shelly remains permits the discri¬ 
mination of rocky and sandy subtidal samples of shelly re¬ 
mains. 

R-mode cluster analysis of species of the living fauna and 
species present in the shelly remains reveals species associa¬ 
tions typical of the major habitats in Bahia la Choya. Analysis 
of species present in the live fauna defines five species associa¬ 
tions: Cerithium — Cardita (rocky outer flats), Tagelus - 


Tbeodoxus (channel and inner flats), Lucina - Oliva (mid to 
outer sand flats), Tellina — A gar on la (shallow subtidal sand 
flats), and Ceritbidea - Mytella (salt marsh). Analysis of spe¬ 
cies in the shelly remains defines 12 clusters, three of which 
are an artifact ol the clustering technique. The remaining nine 
species clusters conform to the area’s major habitats, with 
more than one species association common within each habi¬ 
tat: Cardita - Caryocorbttla and Petricola - Crucibulum 
(both mixed rock and sand flats), Agaronia - Glycymeris, 
Chione californiemis - Felanidla, and Oliva - Parvilucma 
(all three occur in sandy, outer Hat habitats), Chione cortezi — 
Modiolus and Cerithium — Tbeodoxus (both occur in the 
channel), Nassarius — Ceritbidea and Tagelus — Diodora 
(both occur in the high intertidal firmgrounds of the marsh 
and inner flats). Transportation affects the composition ol the 
Tagelus - Diodora cluster. Time-averaging results in high di¬ 
versity and increases the resolution of the cluster analysis of 
the shelly remains. 


KURZFASSUNG 


Die Q-mode Cluster Analyse von Proben der Lebendfauna 
und der Taphozonosen ergibt Gruppierungen, die den wich- 
tigsten Habitaten von Bahia la Choya entsprechen. Die Ana¬ 
lyse von Proben der Lebendfauna definiert vier Gruppierun¬ 
gen: (1) Innenwatt und Salzmarsch, (2) aufteres und mittleres 
Sandwatt, (3) Felswatt und Gezeitenrinne und (4) Subtidal. 
Die Analyse von Proben der Taphozonosen definiert fiinf 
Gruppierungen: (1) Innenwatt, Salzmarsch und Gezeiten¬ 
rinne, (2) sandiges Subtidal, (3) aufieres und mittleres Sand- 
watt, (4) felsiges Subtidal und (3) Felswatt. Postmortaler 
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Transport ist die Ursache fur die Gruppierung von Taphozo¬ 
nosen der Gezeitenrinne mit denen der Salzmarsch und des 
lnnenwatts. Eine langiristige zeitliche Mitteilung (time-aver¬ 
aging) und die Anhaufung diagnostischer Arten in den arten- 
reicheren Proben der Taphozonosen erlauben die Unter- 
scheidung von Proben des felsigen und sandigen Subtidals. 

Die R-mode Cluster Analyse von Arten der Lebendfauna 
und der Taphozonosen ergibt Assoziationen, die die wichtig- 
sten Habitate von Bahia la Choya charakterisieren. Die Ana¬ 
lyse der lebenden Arten fiihrt zu fiinf Assoziationen: Ceri¬ 
thium — Cardita (felsiges AuBenwatt), Tagelus — Tbeodoxus 
(Gezeitenrinne und Innenwatt), Lucina - Oliva (mittleres — 
aufteres Sandwatt), Tellina — Agaronia (flaches sandiges Sub¬ 
tidal) und Ceritbidea - Mytella (Salzmarsch). Die Analyse 
der Arten der Taphozonosen fiihrt zu 12 Gruppierungen, von 
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denen drei ein Kunstprodukt der CIuster-Analyse sind. Die 
ubrigen neun Artengruppierungen entsprechen den wichtig- 
sten Habitaten von Bahia la Chova, wobei jedes Habitat 
durch mehr als eine Assoziation vertreten ist: Cardita — Ca- 
ryocorbula und Petricola — Crucibulum (in Bereichen, in de¬ 
nen Fels- und Sandwatt sich kleinraumig ablosen); Agaronia 

— Glycymeris, Chione californiensis — Felaniella und Oliva 

— Parvilucina (alle drei treten im sandigen Aufienwatt auf); 


Chione cortezi - Modiolus und Ceritbium - Theodoxus (fin- 
den sich in der Gezeitenrinne); Nassarius — Cerithidea und 
Tagelus — Diodora (kommen in Festgriinden des Innenwatts 
und der randlichen Salzmarsch vor). Schalentransport beein- 
fluftt die Zusammensetzungder Tagelus — Diodora Gruppie- 
rung. Zeitliche Mittelung der Taphozonosen fiihrt zu einer 
hohen Diversitat und erhoht die Auflosungsscharfe der CIu- 
ster-Analyse. 


INTRODUCTION 


In this paper we use cluster analysis to define groups of 
samples and species in Bahia la Choya. Cluster analysis is a 
family of techniques used to objectively place entities into a 
hierarchical classification. In paleoecologv, cluster analysis is 
often employed in an elfort to reveal ’’natural” groupings of 
samples or species. When samples are grouped, the geogra¬ 
phic distribution ol the groups may reveal the geographic dis¬ 
tribution of communities. Species groupings help reveal the 
composition of the communities. Davis (1986) and Ander 
berg (1973) provide a good introduction to the mathematics 
of cluster analysis and similar techniques. Cluster analysis 
programs are frequently available within many popular sta¬ 


tistical software packages written for microcomputers and 
mainframes. 

Cluster analysis has been widely employed in taphonomic 
studies. Particularly good examples can be found in Scott & 
West (1976). The technique is often used in order to compare 
the distribution and composition of the live communities to 
the distribution and composition of communities defined on 
the basis of only shelly remains. We continue that tradition in 
this paper. In addition, we illustrate how the taphonomic 
process of time-averaging may actually enhance the resolu¬ 
tion of dendrograms produced by cluster analysis of shelly re¬ 
mains. 


METHODS 


Our analyses are based on samples of the live fauna and 
samples containing shelly remains. A full description of field 
sampling, sample processing, and sample locations is given by 
Fursich et al. (this volume). We considered only those species 
that made up one percent or more of the total number of indi¬ 
viduals in the sample. The data are given in Appendices 1 — 8. 

Forty species were found alive in one or more of 92 sam¬ 
ples. 119 species occurred as shelly remains in one or more of 
101 samples. 

We used the SPSS-X package of programs available on the 
University of Arizona mainframe computer (see Norusis 
1985). 

We used the cosine 0 measure of similarity to calculate 
similarity between each pair of samples or species. The cosine 
0 measure of similarity recalculates abundances as percenta¬ 
ges of the total number of individuals in the sample. 

The principle of the cosine 0 similarity measure is illustra¬ 
ted in Fig. 1. Consider a simple case, where two samples, each 
containing two species are to be compared. The relative abun¬ 
dance of the two species in a sample defines a vector in the 
two-dimensional space illustrated. The degree of similarity 
between these two samples can be measured by the cosine of 
the angle 0 that ist between them. In practice, the number of 
dimensions is equal to the number of species in both samples, 
and cosine 0 is an angle measured in n-dimensional space. 

The similarity measure takes into account both diflerences 
in species composition and differences in the relative abun¬ 
dance of species in common. 



Fig. 1. Example of the cosine 0 similarity index. If Variable 1 is 
the relative abundance of species I and Variable 2 is the relative abun¬ 
dance of species 2, Samples A and B can be considered vectors whose 
position reflects the relative abundance of the species in the two samp¬ 
les. The angle 0 is a measure of the similitarity of the samples. When 
the angle between samples equals 0° (when the samples are identical in 
species composition and relative abundance) then the cosine of 0 
equals 1. When the angle between samples equals 90° (when the samp¬ 
les have no species in common) then the cosine of 0 equals 0. Figure 
from Meldahl (1990). 







We used the conventional unweighted pair-group method 
to form the hierarchical clusters. This methods links samples 
(or species, or groups of samples or species) according to the 
greatest average similarity between samples (or species or 
groups of samples or species). The SPSS-X program rescales 
the similarity measure to a distance measure that varies from 
0 to 25 (labelled "rescaled distance” in Figs. 2 and 3). A small 
distance value indicates high similarity. A high distance value 
indicates low similarity. 
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We use cluster analysis as an exploratory technique. No 
widely accepted method exists to assign statistical significance 
to the resulting clusters. The validity of the resulting clusters 
must be assessed in light of both the appropriateness of the 
methods and the reasonableness of the ecological and paleo- 
ecological interpretations. 


ANALYSIS OF SIMILARITY AMONG SAMPLES 


The cluster analysis of samples groups sample localities 
based on their compositional (e. g. faunal) similarity. This is 
frequently termed a Q-mode analysis. In this mode of analy¬ 
sis, samples of near-identical faunal composition and relative 
abundance w'ill belong to the same group or cluster. Samples 
with no species in common will belong to different clusters. 

We performed two Q-mode cluster analyses: one grouped 
samples based on the similarity of the live fauna; the other 
grouped samples based on the similarity of the shelly remains. 

Q-modeanalysisofsamplescontaininglivefauna 

The dendrogram of samples containing live fauna is shown 
in Fig. 2a. 92 samples contained living species and the clusters 
distinguish four groups of samples. 

1. Innerflats and marsh. These samples come from the 
firmgrounds of the inner flats and the elevated, firm, vege¬ 
tated substrates of the marsh. 

2. A large group of samples clusters at a high level of simi¬ 
larity and groups samples taken in the outer and middle 
flats. These samples are from the sandy portions of the mid 
and low intertidal zone. 

3. The rocky flats and channel cluster groups samples 
taken from areas of hard substrates: either from the rocky 
flats in the northern part of Bahia la Chova or from the shell 
pavements in the main tidal channel in the south. 

4. The sub tidal cluster includes samples taken from the 
shallow subtidal. 

Q-mode analysis of samples containing shelly 
remains. 

The dendrogram of samples containing live fauna is shown 
in Fig. 2b. 101 samples contained shelly remains and the clus¬ 
ters distinguish five groups of samples. 

1. Samples from the inner flats, marsh, and channel 
form a distinct cluster. These samples are from areas of firm 
(inner flats and marsh) or hard substrates and are often domi¬ 
nated numerically by individuals of either Cerithium stercus- 
muscarum or Cerithidea mazatlamca. 

2. Sandy subtidal samples group together. These 
samples come from the shallow subtidal in the sandy, south¬ 
ern of Bahia la Choya. 


3. Samples from the outer and middle flats form a large 
cluster. These samples range from the mid to low intertidal 
and are typically from sandy substrates. 

4. A group of samples from the rocky subtidal forms a 
distinct cluster. These samples are from the rocky, northern 
part of Bahia la Choya. 

5. Samples from the rocky flats comprise the fifth clu¬ 
ster. These samples typically contain encrusting, nestling or 
boring species. Samples come from the hard substrates of the 
northern part of Bahia la Choya and from a few isolated out¬ 
crops in the southern portion of the area. 

Discussion 

Both Q-mode cluster analyses display a common pattern: 
adjacent samples tend to cluster together. This pattern is a 
consequence of the fact that adjacent samples tend to have a 
similar faunal composition. In the cluster analysis of samples 
containing live fauna, this results from similar habitat prefer¬ 
ences. In the analysis of the shelly remains, both original habi¬ 
tat preferences and common response to taphonomic pro¬ 
cesses cause the close association of nearby samples. 

Both analyses reveal habitat-based clusters. Clusters are 
composed of samples from similar substrates and similar du¬ 
rations of tidal exposures. 

The cluster analyses differ somewhat in the specific habitat 
groupings that are formed. Both distinguish a large group of 
samples from the outer and middle flats, an area of rippled 
sand extending from the mid-intertidal to the low intertidal. 
Both analyses serve to distinguish samples from the subtidal 
zone. These samples form a single cluster when based on the 
living fauna. Analysis of the shelly remains serves to distin¬ 
guish samples from the rocky subtidal from those of the 
sandy subtidal. 

Samples from the tidal channel are included with the rocky 
flats samples in the analysis of the live fauna. Analysis of the 
shelly remains, however, places tidal channel samples to¬ 
gether with samples from the inner flats and marsh. The inclu¬ 
sion of channel samples with inner flat and marsh samples in 
the dendrograph of shelly remains may be a consequence of 
the “contamination” of tidal channel samples with many indi¬ 
viduals of Cerithidea mazatlamca that are transported from 
within the marsh. Among samples containing living speci- 
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RESCALED DISTANCE 
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Fig. 2. a. Dendrogram of Q-mode cluster analysis of samples con¬ 
taining live fauna, b. Dendrogram of Q-mode cluster analysis of 
samples containing shelly remains. See text for explanation and dis¬ 
cussion. Figure from FOrsich & Flessa (1987). 
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mens, the rocky flats and channel are characterized by rela¬ 
tively large populations of both Cerithium stercusmuscarum 
and Theodoxus luteofasciatus - grazing gastropods that feed 
on the algae growing on the rocky substrates of the northern 
flats or on the shell pavement that lines the tidal channel. 

Dendrograms from Q-mode cluster analyses of shelly re¬ 
mains are often evaluated according to how well they match 
the dendrographs of analyses based on the living fauna. By 
this qualitative standard, they match reasonably well, dif¬ 
fering principally in the inclusion of channel samples with 
rocky flat samples in one analysis (of the live fauna), and with 
inner flat and marsh samples in the other analysis (of the 
shelly remains). This particular mis-match is evidence of ta- 
phonomic processes. We interpret the clustering pattern in 
the analysis of shelly remains to reflect the transportation of 
Cerithidea mazatlanica shells from the marsh to the channel 
and inner flats. 

In general, the good match is evidence that post-mortem 
processes do not completely redistribute shelly remains in a 
homogeneous fashion. 


The cluster analysis of shelly remains may also be evaluated 
according to its power of resolution. By this standard, the 
cluster analysis does quite well because it serves to distinguish 
two groups of subtidal samples, rather than the one revealed 
in the analysis of the live remains. The greater resolution of 
habitats possible with shelly remains is a consequence of the 
much higher diversity of shelly remains. The greater number 
of species present among the shelly remains results from ex¬ 
tensive time-averaging. Rare or transient species that may be 
diagnostic of particular habitats are more likely to be found in 
time-averaged accumulations of shells. 

Our interpretation and comparison of the two dendro¬ 
grams reveals that one taphonomic process (transportation) 
may serve to obscure important habitat distinctions — in the 
case of the channel samples. However, another taphonomic 
process (time-averaging) serves to enhance the resolution of 
habitats - in the case of the discrimination of the subtidal 
samples. This latter phenomenon is a good example of how 
taphonomic effects may add important information to a 
paleoenvironmental analysis. 


ANALYSIS OF SIMILARITY AMONG SPECIES 
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The cluster analysis of species occurrences groups species 
together based on the degree to which they co-occur in sam¬ 
ples. This is frequently termed an R-mode analysis. In this 
mode of analysis, species that invariably occur in the same 
samples will belong to the same group or cluster. Species 
which never occur in the same sample will be found in dif¬ 
ferent clusters. 

We performed two R-mode cluster analyses: one grouped 
living species based on their co-occurrences among the sam¬ 
ples; the other grouped species of shell remains based on their 
co-occurrences among our samples. 

R-mode analaysis of species found alive in 
samples 

The dendrogram of the 40 species found alive in one or 
more of our samples in shown in Fig. 3. The clusters form at 
relatively low levels of similarity, forming five species asso¬ 
ciations. 

1. The Ccrithium - Cardita association is characterized by 
epifaunal gastropods and nestling or encrusting bivalves. 
Members of this association are most common in the mixed 
rock and sand flats of the low intertidal in the northern part of 
the bay. 

2. The Tagelus - Tbeodoxus association also includes two 
common shallow-infaunal venerid bivalves, Chione califor- 
niensis and Protothaca grata and the epifaunal gastropod Nas- 
sarius brunneostoma. These species typically occur in the vici¬ 
nity of the main tidal channel and in the firmgrounds of the 
middle to inner flats. 

3. The Lucina — Oliva association is dominated by the 
small infaunal bivalves Fclaniella, Lucina lampra , Lucina 
rnazatlanica, and Parvilucina. The gastropods Oliva spicata 
and Nassarius iodes are also typical. This assemblage is typical 
of the mid to outer sandflats of the low intertidal. 

4. The Tellina - Agaronia testacea association is charac¬ 
terized by the shallow-burrowing echinoid Encopc micropora 
and the shallow-burrowing gastropod Agaronia testacea. 
Members of this association typically occur in the shallow 
subtidal and lowest intertidal of the outer flats. 

5. The Cerithidea — Mytella association consists of three 
species typical of the marsh/estero. 

R-mode analysis of species found as shelly 
remains in samples 

The dendrogram of the 119 species present in one or more 
sample of the shelly remains is shown in Fig. 4. The clusters 
form at relatively high levels of similarity, forming twelve 
species associations. 

1. The Tellina — Pteria association is dominated by shal¬ 
low burrowing and epifaunal bivalves. Many of the species 
are characteristic of the shallow subtidal of the mixed sandy 
and rocky flats of the northern part of Bahia la Choya. 

2. The Cardita — Caryocorbula association is composed of 
dominantly hard substrate bivalves (either nestlers or cemen- 


ters) and epifaunal gastropods. This association is often found 
in the low intertidal and shallow subtidal mixed rock and sand 
flats of the northern part of Bahia la Choya. 

3. The Petricola — Crucibulum association is also com¬ 
posed of bivalves and gastropods typical of hard substrates 
(e. g. the two name-giving genera), plus a few species of small, 
shallow-burrowing bivalves. This association occurs in the 
mixed rock and sand flats in the shallow subtidal of the north¬ 
ern part of Bahia la Choya. 

4. The Chione cortezi — Modiolus association is most 
frequently found in along the main tidal channel. 

5. The Chione californiensis — Felaniella association is 
dominated by shallow-burrowing bivalves that typically 
occur in the mid to outer sand flats of the southern portion of 
Bahia la Choya. 

6. The Oliva — Parvilucina association is dominated by 
species typical of the middle to outer sand flats. The inclusion 
of the hard substrate species Vermetus and Crassostrea, how¬ 
ever, is a puzzle. 

6 a. Nassarius brunneostoma and Cerithidea rnazatlanica 
form a distinct subgroup with the Oliva - Parvilucina clus¬ 
ter. N. brunneostoma and C. rnazatlanica typically occur in 
inner flat areas. 

7. The Cerithium - Theodoxus association. These two 
gastropods, together with Nassarius iodes , are commonly 
found in the main tidal channel and on the firmgrounds of the 
inner flat. 

8. The Agaronia — Glycymeris association consists of spe¬ 
cies typical of the outer sand flats. 

9. The Anomia — Crepidula association is dominated by 
encrusting species. This assemblage is suspect because it is 
defined principally by co-occurrence within a single subtidal 
sample. 

10. Species in the Dosima — Chama association are typi¬ 
cally found in the lowest intertidal and shallow subtidal of the 
sandv and rocky outer flats. However, this association is 
defined largely by the influence of a single, subtidal sample. 

11. The Tagelus — Diodora association is a mixture of spe¬ 
cies inhabiting the inner flats near vegetated areas ( Tagelus 
and Gregariella ) and those washed out of the spit deposits at 
the strand line ( Diodora and Agladrillia). 

12. The Isognomon — Crassispira association is an outer- 
flat, shallow subtidal assemblage found on the mixed rock and 
sand substrates of the northern part of Bahia la Choya. This 
assemblage is suspect because it is defined by the co-occur¬ 
rence of three species in a single sample. 

Discussion 

The most dramatic difference between the two dendro- 
graphs is their size. Only 40 species were found alive in suffi¬ 
cient abundance in one or more sample. Almost three times as 
many (119) occurred in sufficient relative abundance (more 
than 1 %) in samples of shelly remains. The difference in di¬ 
versity is a consequence of time-averaging. The very low sedi¬ 
mentation rates in Bahia la Choya permit the accumulation of 
the shells of rare and infrequent species. 
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The greater number of species among the shelly remains 
also includes several species found only in one sample. These 
species are those that do not form subgroups and join the den- 
drograph only at very low levels of similarity. These are the 
species listed at the top of the dendrogram. 

When two or more species co-occur in a single sample but 
occur in no other samples, the resulting high similarity among 
the species may result in a distinct cluster. These single¬ 
sample clusters (e. g. the Dosinia — Chama and Isognomon — 
Crassispira associations described above) should be consid¬ 
ered suspect. Such clusters are based on rare species that occur 
in only one sample and should not be accorded much signifi¬ 
cance. In a sense, they represent accidents of species distribu¬ 
tions and sampling. 

The greater number of associations distinguished by the 
analyses of the shelly remains is a simple consequence of the 
greater diversity of species found among the shelly remains. 
More species (119 vs. 40) in more samples (92 vs. 101) will 
tend to produce more clusters. 

Species in the living fauna form five clusters which roughly 
correspond to major habitats within Bahia la Choya. The Ce- 
ritbium — Cardita association is typical of the hard substrates 
of the outer flats and the shell pavements of the tidal channel. 
The Tagelus — Theodoxus association is common along the 
tidal channel and inner flats. The middle and outer sand flats 
fauna typically includes members of the Lucina — Oliva spi- 
cata association. The Tellina — Agaronia association occurs 
in the shallow subtidal. The Cerithidea — Mytella association 
characterizes marsh habitats. 

Associations formed by species found among the shelly re¬ 
mains also conform to Bahia la Choya’s principal habitats. In 
many cases, however, more than one association occurs with¬ 
in a single habitat. Sandy, outer flat habitats contain three 
associations ( Agaronia — Glycymeris , Cbione californiensis — 


Felamella and Oliva — Parvilucina ), and the mixed rock and 
sand habitat of the low intertidal and subtidal includes three 
associations ( Tellina — Pteria , Cardita — Caryocorbida and 
Petricola — Crucibulum). Two associations are commonly 
found in the channel ( Chione cortezi — Modiolus and Ceri- 
tbium — Tbeodoxus). The high intertidal firmgrounds of the 
inner flats include two associations (A ’assarius — Cerithidea 
and Tagelns — Diodora ). We dismiss the environmental signi¬ 
ficance of the remaining associations because they are often 
based on a single sample. 

We interpret the greater resolution of the dendrogram of 
shelly remains as a result of the time-averaged nature of the 
shelly deposits. The long time during which shells have accu¬ 
mulated has allowed the accumulation of species diagnostic of 
the area’s habitats and sub-habitats. Many of these species 
were probably not contemporaneous, but were nevertheless 
characteristic of a relatively limited set of environmental con¬ 
ditions. Samples of species that happen to be alive during a 
particular sampling interval are not likely to include all of the 
species that might inhabit the habitat over an extensive time 
interval. We are tempted to argue, therefore, that — as long as 
no environmental changes are involved — time-averaged 
samples provide a better means by which to discriminate 
habitats than samples that represent an “instant” in time. This 
would be the case in general only if between-habitat transpor¬ 
tation was negligible. 

The effects of transportation are evident in some of the as¬ 
sociations formed from species in the shelly remains. For ex¬ 
ample, the Tagelus — Diodora association contains species 
that are winnowed out at the strandline of the spit. Thus, spe¬ 
cies from the spit deposits are recycled back into active sedi¬ 
mentary environments. The Oliva — Parvilucina association 
includes some species more characteristic of rocky substrates. 
Those species may have been transported in form the adjacent 
rocky flats. 


CONCLUSION 


We conclude that cluster analysis of shelly remains is an 
effective exploratory tool in paleoenvironmental analysis. 
A good correspondence exists between clusters based on the 
distribution of live fauna and the clusters based on the distri¬ 
bution of shelly remains. 

The composition and distribution of the Bahia la Choya 
taphocoenoses is a function of both original species distribu¬ 
tions, the original species composition of habitats and the 
taphonomic processes of time-averaging and transportation. 
Of the two taphonomic processes, time-averaging appears to 
exert the greatest effect. One beneficial side effect of time¬ 
averaging is the greater resolution achieved with cluster ana¬ 


lysis of shelly remains. Time-averaged samples of shelly re¬ 
mains are more likely to contain environmentally diagnostic 
species. 

Because both species composition and taphonomic proces¬ 
ses vary significantly from habitat to habitat, the geographic 
distribution of the taphocoenoses, as shown in the dendro¬ 
gram, is an excellent indicator of important habitat differen¬ 
ces within Bahia la Choya. With the exception of species asso¬ 
ciations resulting from co-occurrence in single samples, spe¬ 
cies associations defined by cluster analysis permit the rec¬ 
ognition of environmentally sensitive taphocoenoses. 


◄ 


Fig. 4. Dendrogram of R-mode cluster analysis of species found as shelly remains in one or more samples. 
Italicized names indicate names given to species associations. See text for explanation and discussion. 



88 


REFERENCES 


ANDERBERG, M. R. (1973): Cluster Analysis for Applications. - 
359 pp., Academic Press, New York. 

Davis, J. C. (1986): Statistics and Data Analysis in Geology. — 
646 pp., John Wiley and Sons, New York. (2nd edit.) 

FORSICH, F. T. & FLESSA, K. W. (1987):Taphonomy of tidal flat mol¬ 
luscs in the northern Gulf of California: Paleoenvironmental 
analysis despite the perils of preservation. - Palaios, 2, 
543-559, Tulsa, OK. 

Fursich, F. T„ Flessa, K. W., Aberhan, M., Feige, A. & SchOdl- 
BAUER, S. (1991): Sedimentary habitats and molluscan faunas of 
Bahia la Choya (Gulf of California, Sonora, Mexico). - Zitte- 
liana, 18, 5 — 51, 10 pis, Munchen. 


MEHLDAHL, K. H. (1990): Paleoenvironmental and stratigraphic im¬ 
plications of taphonomic processes: Case studies from Recent 
and Pleistocene marine environments. — 217 pp., unpublished 
Ph. D. Dissertation, University of Arizona, Tucson, AZ. 

NORUSIS, M. J. (1985): SPSS-X Advanced Statistics Guide. - 505 pp., 
McGraw-Hill, New York. 

Scon, R. W. & West, R. R. (eds) (1976): Structure and Classification 
of Paleocommunities. — 291 pp., Dowden, Hutchinson and 
Ross, Stroudsburg, Pennsylvania. 


